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ABSTRACT 
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5^ ; Forty molecular cloud cores in the southern hemisphere from the initial Spitzer 

Space Telescope Cores-to-Disks (c2d) Legacy program source list have been sur- 
veyed in 13 CO(2^1), 12 CO(4^3), and 12 CO(7^6) with the Antarctic Submil- 
limeter Telescope and Remote Observatory (AST/RO). The cores, ten of which 
contain embedded sources, are located mostly in the Vela, Ophiuchus, Lupus, 
\Q • Chamaeleon, Musca, and Scorpius complexes. 12 CO(7— >6) emission was un- 

detected in all 40 clouds. We present data of 40 sources in 13 CO(2^1) and 
' 12 CO(4^3), significant upper limits of 12 CO(7— >6), as well as a statistical anal- 

ysis of the observed properties of the clouds. We find the typical 13 CO(2— >1) 
linewidth to be 2.0 km s -1 for cores with embedded stars, and 1.8 km s _1 for 
all others. The typical 12 CO(4^3) linewidth is 2.6 to 3.7 km s -1 for cores with 
known embedded sources, and 1.6 to 2.3 km s _1 for all others. The average 13 CO 



column density derived from the line intensities was found to be 1.9 x 10 15 cm -2 
for cores with embedded stars, and 1.5 x 10 15 cm -2 for all others. The average 
kinetic temperature in the molecular cores, determined through a Large Velocity 
Gradient analysis of a set of nine cores, has an average lower limit of 16 K and 
an average upper limit of 26 K. The average molecular hydrogen density has an 
average lower limit of 10 2 ' 9 cm -3 and an average upper limit of 10 3 ' 3 cm -3 for 
all cores. For a different subset of nine cores, we have derived masses. They 
range from 4 to 255 M . Overall, our c2d sample of southern molecular cores 
has a range of properties (linewidth, column density, size, mass, embedded stars) 
similar to those of past studies. 
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1. INTRODUCTION 

Star forma tion occurs in dense molecular clouds and their dense co res are the progenitors of 
protostars (jdi Francesco et al. 1120061 ; IWard- Thompson et al. II2006I ). Such cores are compact 
areas (< 0.5 pc) of molecular gas at relatively high density (> 10 3 cm -3 ). Large efforts have 
been made in the observation and characterization of such cores. One of the biggest and 
most recent of these efforts is the Spitze r Space Telescope ( SST) program 'From Molecular 
Cores to Planet-forming Disks' or 'c2d' (jEvans et al. II2003I ). The c2d program surveyed at 
infrared wavelengths over 100 relatively isolated cores that span the evolutionary sequence 
from starless cores to protoplanetary disks (Huard et al. 2006). These sources include a 
wide range of cloud masses and star-forming environments. In order to understand how 
cores are formed and subsequently evolve into stars, it is necessary to understand the gas 
dynamics of the underlying molecular cloud. The primary means of studying molecular gas 
are observations of emission by the rotational transitions of the ground vibrational state 
of carbon monoxide (CO). These spectral lines occur at frequencies of J x 115 GHz, for 
transitions from the J to the J — 1 state. A realistic picture of the thermodynamic state of 
a molecular gas can be gained through observations of several transitions of CO isotopes. If 
the brightness of several of these lines is known, models of the r adiative transfer, such as a 
Large Velocity Gradient (LVG) model (IGoldreich fc Kwan Ill974 ). can be used to determine 
density and temperature of the gas. In order to avoid degeneracy in this solution, it is 
advisable to have observations of a high- J transition line that is sufficiently high in energy 
to be only weakly populated. 

From the initial c2d source list of Evans et al. (2003), we have observed with AST/RO 
all 40 molecular cores that are visible from the South Pole (5 < —16°), in the 13 CO(2-^l), 
12 CO(4— >3), and 12 CO(7— >6) transitions. (The initial SST c2d source list was later shortened 
due to the time available to the c2d program. Some of the sources in this paper have not 
been observed with the SS T.) Ten out of the 40 cores are known to have sources embedded in 
them according to IRAS (jEvans et al. 1120031 ; iLee fc Myers Ill999l ). The detect ion threshold 
for st ar formation using IRAS is typically taken to be ~0.1 L (d/140pc) (IMyers et al. 
19871 ). Cores in which IRAS has not detecte d a source, which in this paper are referred to 



as 'starless', may s t ill contain faint sourc es (jYoung et al. I l2004j ; iBourke et al. 1 120051 . 12006 



Huard et al. 1 120061 ; iDunham et al. II2006I ). The c2d team is producing a list of candidate 



very low luminosity objects (VeLLOs) in all cores. The multiple transitions observed allow us 
to perform an LVG analysis and gain information about the kinematic and thermodynamic 
conditions in a large sample of molecular cores. These observations and results can be used 
in conjunction with the c2d infrared results for further studies of these cores. 



In §2, we describe the observatory, the observations, and the method of data reduction. §3 
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presents the dataset in the form of spectra, maps, tables, and a statistical analysis. The 
Large Velocity Gradient analysis of the radiative transfer is presented in §3.1, which allows 
a determination of kinetic temperature and molecular hydrogen density. Masses for several 
cores are derived in §3.2. In §4, we summarize our results and present our conclusions. 



OBSERVATIONS 



The observations were performed during the austral winter of 2005 at the Antarctic Sub- 
millimeter Telescope and Remote Observatory (AST/RO), which was located at 2847 m 
elevation at the Amundsen-Scott South Pole Station. These are among the final obser- 
vations of AST/RO prior to its decomissioning in December 2005. Low water vapor and 
high atmospheric stability make this site exceptionally good for submillimeter astronomy 



( jChamberlin et al. 1 119971 ; lLane I Il998l ). AST/RO is a 1.7 m diamete r, offset Gregorian 



telescope that obse rves at wavelengths between 1.3 mm and 200 /im (IStark et al. 1 12001 



Oberst et al. 1120061 ). Emission from the 220.399 GHz J = 2 -> 1 transition of 13 CO 
was map ped using an SIS receiver with 80 to 160 K double-sideband receiver noise tern- 



perature ( 


Kooi et al. 


1992) 


Honineh et al. 


1997) 


was i 



A dual-channel SIS waveguide receiver ( IWalker et al. 1 11992 
19971 ) was used for simultaneous observation of 461 GHz and 807 GHz. 
Single-sideband receiver noise temperatures were 340 to 440 K during observations of the 
461.041 GHz J = 4 -> 3 transition of 12 CO and 700 to 1100 K for the J = 7 -> 6 transition 
of 12 CO at 806.652 GHz. 

A multiple-position-switching mode was used with emission-free reference positions at least 
30' from the mapped area. Integration times per point were 347 seconds at 220 GHz and 
220 seconds at 460/807 GHz. The beamsizes (FWHM) were 180" at 220 GHz, 103" at 460 
GHz, and 58" at 807 GHz. All maps were fully-sampled with neighbouring points 1.5' apart 
at 220 GHz and 0.5' apart at 461/807 GHz. Tests of the AST/RO pointing model indicate 
that pointing errors at the time these data were taken could in some cases be as large as 2' 
at 220 GHz and V at 461/807 GHz. 

The telescope efficiency r]e was estimated through the reduction of a series of skydips to ~ 
81% for the 220 GHz receiver, « 72% at 461 GHz, and « 80% at 807 GHz. As described in 
detail in Stark et al. (2001), the AST/RO main beam efficiency and telescope efficiency are 
identical to the level that we can measure. The telescope efficiency has been applied to the 
data as presented below. Atmosphere-corrected system temperatures at the three frequencies 
ranged from 280 to 480 K, 100 to 2000 K, and 4000 to 10000 K, respectively. Two acousto- 
optical spectrometers (AOS) (jSchieder et al. Ill989l ) with 670 kHz channel width, 1.1 MHz 
resolution bandwidth, and 1.5 GHz total bandwidth were used as a backend, resulting in 
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velocity resolutions of 0.91 km s 1 at 220 GHz, 0.65 km s 1 at 461 GHz, and 0.37 km s 1 
at 807 GHz. For obse rvationos at 13 CO(2^1) we used a higher resolution AOS (HRAOS) 



(ISchieder et al. 1119891 ) in addition to the two mentioned above. The HRAOS has a channel 
width of 31.6 kHz and a resolution bandwidth of 64 MHz, resulting in a velocity resolution 
of 0.043 km s _1 at 220 GHz. The data from this spectrometer are used to derive linewidths 
when the lower resolution was not sufficient and to confirm linewidths close to the resolution 
limit at 220 GHz. These data are also marked in the tables. 

A chopper wheel cal ibration technique was used to observe the sky and a blackbody load of 



known temperature (IStark et al. 1120011 ) every few minutes during observations. Skydips were 
performed frequently to determine the atmospheric transmission. To monitor the pointing 
accuracy, bright sources (NGC3576 and/or the moon) were mapped every few hours. After 
regular intervals or after changes to the receiver systems (tuning), the receiver was manually 
calibrated against a liquid nitrogen load and a blackbody load of known temperature. This 
process also corrects for the dark current of the AOS's optical CCDs. The internal consistency 
of the AST/RO flux scale was estimated from a large-scale map of Lupus, also taken in 2005 
with identical instrumental settings. This large map was constructed from smaller maps, 
with some overlap between them, yielding over a thousand points towards which more than 
one spectrum was measured at different times. From these data, the flux scale is reproducable 
at the 30% level after correction for atmospheric absorption (Tothill et al., in preparation). 

The data taken for this survey were frequency calibrated in several ways. Before and during 
observations, we observed a comb spectrum with every automatic receiver calibration to 
monitor the frequency scale over the entire observing period. During measurements of the sky 
performed with every automated calibration, the telescope observes a mesospheric telluric 
line. The velocity scale of all observations was referenced to this scale, resulting in an 
accuracy of plus/minus one channel. 



he d ata were reduced using the COMB data reduction package as described in lStark et al 



(120011 ). Linear baselines were removed in all spectra at all frequencies, excluding regions 
where emission was expected base d on previous 13 CO (J = 1 — > 0) observations (R. Otrupcek, 



private communication; see also lOtrupcek et al.l (120001 )). Removal of first-order baselines 



proved to be generally sufficient for the 13 CO(2^1) and 12 CO(4— >3) data. In very rare cases 
(0.01% of all data), two or three Fourier components had to be removed in addition to the 
linear baseline. At 12 CO(7— >6) we did not detect significant (> 3d) emission. However, we 
were able to deduce upper limits for the 12 CO(7— >-6) intensity from our data for all of the 
observed sources. Since we detected strong 12 CO(7— >6) emission in our pointing sources, 
we can be certain that the non-detection of this line in the cores is not a result of receiver 
failure. We conclude that the 12 CO(7^6) emission from the observed cores is so faint that it 



- 5 - 



lies below the detection sensitivity of our detectors. We thus deduce significant upper limits 
for this emission. 

A list of observed sources in order of increasing R.A., coordinates of the map centers, as well 
as mapsizes is given in Table [TJ Cores with known embedded IRAS sources are marked (*). 
The molecular cores w e observed contain all cores from the source list of the Spitzer Space 
Telescope c2d project (jEvans et al. 1120031 ) that are visible from the South Pole (5 < —16°). 
The observed mapsizes were chosen to cover at least the areas observed by c2d, but are 
larger in almost all cases. 12 CO(4^3) mapping was sometimes confined to regions where 
previously mapped 13 CO(2-^l) emission was observed. 



3. RESULTS AND ANALYSIS 

For each observed core, we present the highest intensity spectrum of both detected transitions 
in Figures la - le. Observed and derived parameters, and source information are listed 
in Table [2] for the 13 CO(2^1) transition and in Table [3] for the 12 CO(4^3) transition. 
Molecular cores with known embedded IRAS sources are marked with a (*). The R.A. and 
Dec. coordinates in Tables [2] and [3] represent the position of the spectrum with the maximum 
integrated intensity in each map. The source velocity vlsr wa s obtained by fitting a Gaussian 
to the line with maximum intensity. Av corresponds to the FWHM of this fit. In cases where 
the linewidth is lower or close to the spectrometer's velocity resolution of 0.91 km s -1 , we 
confirmed the linewidth from our higher resolution data (0.043 km s" 1 ). Cores for which we 
did this are marked in the table. J Tmb{v) dv gives the intensity integrated over the brightest 
line in the map, and rms is the root-mean-squared noise of the channels in the spectrum. 

In Table [21 the 13 CO(2^1) intensity is converted to a total column density of 13 CO (Ni3co) 
by assuming local thermodynamic equilibrium (LTE): 

M i si v mi4 T ** e¥ Jr(v)dv 

^Vi3 CO = 1.51 X 10 (1J 

I — e r ex 

Assuming optically thin 13 CO (2— >1) emission, that the source fills the main beam, that the 
background temperature of 2.7 K can be neglected and that the Rayleigh- Jeans correction 
can be neglected as well, we use the relation 

T ex t(v) = T MB (v) (2) 

where t(v) is the optical depth of the 13 CO(2— >1) line, Tmb is the main beam brightness 
temperature and T ex is the excitation temperature, which is usually taken to be the kinetic 
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temperature in the cloud: T ex =Tki n =W K (Up to T ex 
on the 13 CO column density is less than a factor of 2). 



40 K the effect of this parameter 



In Figure [2J we show histograms of the linewidths of the two detected transitions. The 
distribution of the 13 CO(2— >1) linewidth has a mean of 1.8 km s _1 and a standard deviation 
of 0.4 km s _1 . The cores with embedded sources are represented by the black bars. It 
is noticeable that their distribution is generally skewed to higher values. The mean of 
their linewidths distribution is 2.0 km s _1 and the standard deviation is 0.4 km s -1 . Cores 
with embedded stars tend to be warmer and more turbulent, and their emission lines are 
therefore subject to more broadening, than are cores without embedded stars. If the gas 
kinetic temperature is typically 20 K, the corresponding optically thin thermal linewidth 
of CO is only 0.2 km s _1 , so the observed widths are due primarily to turbulent motions 
and optical depth effects. DC253.8-10.9 shows an exceptionally narrow line with only 0.8 
km s _1 . It also is the weakest core in 13 CO(2— L328 has the broadest line of the cores 
with no embedded IRAS source at 2.7 km s -1 . This core has recently been found to have a 
VeLLO embedded (Lee, C.-W., in preparation). T he widest line at 2.9 km s _1 is emitted by 
DC275.9+1.9, which does have an embedded star (iBourke et al. Ill995l ). 

The distribution of the 12 CO (4—^3) linewidth also shows that the lines from cores with 
embedded sources tend to be broader than lines from starless cores. The mean linewidth of 
cores with embedded stars is 3.2 km s -1 , with a standard deviation of 1.3 km s _1 . Cores 
with no embedded stars have a mean linewidth of 2.0 km s" 1 , with a standard deviation of 
0.7 km s -1 . For this histogram, we did not consider cores whose peak intensity is lower that 
5a. These omitted cores are DC302.6-15.9, L63, CB68, and L328. The wi dest line of al l 
cores is found in DC303.8-14.2 at 5.8 km s" 1 , which has an embedded star ( Bourke et al. 



19951 ). The widest line of cores with no embedded sources is emitted by DC300. 7-1.0 with 
3.5 km s" 1 . 



The survey of southern cores by lVilas-Boas et al. I (119941 ) in 13 CO(l— > 0) reports an average 
linewidth of 0.8 km s~ l . The only source obs erved in both the Vilas-Boas et al. survey and 
this survey is CG30. IVilas-Boas et al. I (I1994J ) gives a linewidth of 1.19 km s -1 obtained with 
a resolution of 0.12 km s _1 . Using our high resolution spectrometer, we observe d a linewidth 
of 1.2 ± 0.1 km s" 1 , which is in good agreement with IVilas-Boas et al. I (I1994J ) . 



A histogram of the 13 CO column density is shown in Figure [31 The mean 13 CO column 
density for cores with no embedded stars is 0.5 to 2.0 x 10 15 cm -2 and slightly higher, 
around 1.0 to 2.0 x 10 15 cm -2 , for cores with embedded stars, which do not show 13 CO 
column densities lower than 1.0 x 10 15 cm" 2 . Fifty-five percent of the cores have a 13 CO 
column density between 0.5 to 1.5 x 10 15 cm -2 , 



and 75% have column densities between 



0.5 to 2.0 x 10 15 cm" 2 . Assuming a 13 CO/H 2 ratio of 1.7 x 10~ 6 flWilson fc Roodlll994 
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Frerking et al. 1119821 ). we can convert the CO column densities to the H2 column density. 



The H 2 column density distribution has a mean of 1.9 x 10 21 cm" 2 and a standard deviation 
of 0.7 x 10 21 cm" 2 for cores with embedded sources. For all other sources the mean is 1.5 
x 10 21 cm -2 and the standard deviation is 0.7 x 10 21 cm -2 . Cores with embedded sources 
tend to have higher column densities. Seve n of the observ ed cores lie in Ophiuchus which 
has a high level of star formation activity (IWilking 1 1 19921 ). In these cores the H2 column 
density assumes high values, with a mean of 3.7 x 10 21 cm -2 and a dispersion of 1.5 x 
10 21 cm -2 . Twenty-three cores lie within clouds where the star-forming activity is low (Vela, 
Musca, Coalsack, Chamaeleon II and III, Lupus, Scorpius). These cores have a mean H2 

This difference 



column density of 2.7 x 10 21 cm 2 with a dispersion of 1.0 x 10 21 cm 2 . 



as an indicator of star-forming activity was also observed by IVilas-Boas et al. I (120001 ) in a 
13 CO(1-0) survey of different ( none of the source s in th eir and our survey overlap) dense 
condensations in dark clouds. IVilas-Boas et al. I (120001 ) set the limit for star-formation- 
enabling H2 column density to as high as 1.1 x 10 22 cm -2 . This limit appears to be too 
high, considering our results for B59. B59 in Ophiuchus is the strongest core observed and 
has the exceptionally high 13 CO column density of 4.3 x 10 15 cm" 2 and H2 column density 
of 7.3 x 10 21 cm" 2 , thus indicating strong star formation activ ity. A group of embedded 
sources has already been found in this core (IBrooke et al. 1 120071 ). 



3.1. Large Velocity Gradient Analysis 

For a more detailed analysis of our data, we chose only the cores most suitable for the 
purpose: cores with strong (> 5cr) emission in both 13 CO(2^1) and 12 CO(4^3). We also 
made sure there was no large pointing offset between intensity peaks at the two wavelengths. 
In Figures 4a to 4e, we present maps of the nine cores chosen for the LVG analysis. The 
AST/RO contour maps are sup erimposed on optical i mages from the Digital Sky Survey 



(DSS) retrieved from Sky View ( McGlynn et al. Ill996l ). The cores are shown in order of 



ascending R.A. We present the majority of cores superimposed on 15' x 15' images, but 
AST/RO map sizes vary. Maps in the same row show the same source, the left column shows 
the 13 CO(2— >1) emission, the right column shows the 12 CO(4— >3). All maps are labeled with 
the source name, the observed transition, and the maximum integrated intensity in K km 
s" 1 . Contours in all maps are in percentage of the peak intensities. We show contour levels 
in steps of 10% of the highest contour, whose value is given in each figure. The first maps on 
each page show the AST/RO effective beamsizes of 216" at 220 GHz and 117" at 461 GHz. 
Since the maps were made using Gaussian interpolation, the effective beamsizes in the maps 
are slightly larger than the actual telescope beamsizes. 
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Maps that are larger than 15' x 15' are shown in order of ascending R.A. after the set of 
15'x 15' maps. The top figure shows the 13 CO(2^1) map, the bottom figure shows the 
12 CO(4— >3) map. In cases with more than one source per map, the sources are labeled 
individually. 

We can estimate the kinetic temperature Tki n , and the number density of molecular hy- 
drogen n(H?), using a La rge Velocity Gradient (LVG) analysis of the radiative transfer 
(IGoldreich fc Kwan 111974 ). Our LVG radiativ e transfer code has been developed and ap- 
plied by M. Yan and S. Kim (]Kim et al. 1120021 ) and has been implemented as a java applet 
q T his mode l simul ates a plane-parallel cloud geometry and uses the CO collisional rates 
from iTurner I (119951 ). newly-derived values for the H 2 ortho-to-para rati o (~ 2), and the 



collisional quenching rate of CO by H 2 impact (iBalakrishnan et al. 1 120021 ). The model has 
two input parameters: the ratio of 12 CO to 13 CO abundance, and the ratio X(CO)/W, 
where X(CO) is the fractional CO abundance and W denotes the veloci ty gradient. The 



abun dance ratio 12 CO/ 13 CO is taken to be 60 in dense molecular cl ouds (IWilson fc Rood 



199J). The 12 CO/H 2 ratio is taken to be 10" 4 (IFrerking et al. Ill982f ) (leading to a 13 CO to 



H 2 ratio of 1.7 x 10 6 ) and the velocity gradient across the cores 1 km s 1 pc 1 . 

The 12 CO(4— >3) data were convolved with a Gaussian kernel to the lower 13 CO(2- 
lution. The resulting peak intensities are listed in Table HI 



►1) reso- 



Although we were unable to detect any significant 12 CO(7^6) emission, we can use these 
data to derive a significant upper limit for this emission for each source. We convolved the 
higher resolution 12 CO(7— >6) data with a Gaussian kernel to the 13 CO(2— >1) resolution and 
derived a 3a upper limit for the emission through this relation: 



3a = 3 x 



rms x v N x channelwidth x 0.94 
FWHM(4 - 3) 



(3) 



where 'rms' is the rms noise of the co-added spectrum. N is the number of channels over which 
the rms was measured (twice the expected linewidth of the 12 CO(7— >6) line). 'Channelwidth' 
is the width of the acousto-optic spectrometer channels at this wavelength. FWHM(4-3) is 
the half-power width of the corresponding 12 CO(4— >3) line. The factor of 0.94 comes from 
the relation for Gaussian profiles: FWHM x peak = 0.94 x area. 

The results are shown in Table HJ where upper and lower limits are denoted u.l. and 1.1. 
Since the line brightness temperature Ty^g of the 12 CO(7^6) line is an upper limit, the 

To 



ratio T^ 6 



/T 4 ^ 3 , the kinetic temperature and the H 2 density are also upper limits. 



^ttp : / / arcsec . se jong . ac . kr/^skim/lvg2 . html 
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obtain the lower temperature and density limits, we assume no 12 CO(7— >6) emission at all: 
T 12 - 

As can be seen from Table HI the lower temperature limits average to 16 K and are not 
greater than 20 K for any source. The distribution of the upper temperature limits have an 
average of 26 K. The H2 density distribution is fairly uniform, ranging from 10 27 cm -3 to 
10 3 ' 7 cm -3 with an average lower limit of 10 29 cm -3 and an average upper limit of 10 3 ' 3 cm -3 
for all sources. Due to the small sample size, we cannot detect a significant difference in the 
mean densities of cores with and without embedded stars. 



3.2. Masses 



Nine cores were chosen for the derivation of their masses. We selected cores whose half-peak 
contours in 13 CO(2— >1) are closed to derive the cores' masses. These cores are also resolved. 
Maps of the selected sources are presented in Figures 4 and El We show contour levels in 
steps of 10% of the highest contour, whose value is given in each figure. DC302.6-15.9 and 
DC303.8-14.2 are at very high elevation (79° and 77°) and possibly suffer from a systematic 
pointing offset as suggested from the background image. For the calculation of their masses, 
a potential pointing offset is irrelevant. There is no evidence that other cores show a similar 
pointing offset. We can obtain the total mass of the core using 



M = {m) / N(H 2 ) dA 

J A 



(4) 



where A is the projected area and (m) is the mean mass per hydrogen molecule. This mean 
mass per hydrogen molecule (m) accounts for a helium abundance of one helium atom per five 
hydrogen molecules. We take into account that the sources are not spherical, but elliptical, 
and list their two radii (major and minor axes) in Table [5) We derive an average 13 CO column 
density for the core area that is defined by the 50% peak intensity contour by co-adding all 
data that were taken in this area and making use of equation (1). This conversion assumes 
an excitation temperature of T ex = 10 K. Up to T ex = 40 K, the effect of this parameter 
on the 13 CO column de nsity is less than a fac t or of 2. By assuming a ratio of 1.7 x 10 -6 
between 13 CO and H2 (IWilson fc Rood Ill994l ; iFrerking et al. Ill982l ). we then convert the 
13 CO column densities to the H 2 column densities listed in Tableland use the relation above 
to determine their masses. Their distances are quoted from the literature: DC259.5-16.4, 
CG30 DC267.4-7.5, DC253.6+ 2.0, and DC274.2-0.4 lie within the Vela complex at 450 ± 
50 pc (|Woermann et al. 112001); DC 297.7-2.8 is associated with the Coalsack at 150 ± 30 pc 



(IFranco 



1995 



Corradi et al. 



19971 ): DC302.6-15.9 an d DC303.8-14.2 lie in t he Chamaeleon 



cloud complexes at distances between 150 - 180 pc (IWhittet et al. 1 119971 ; iKnude &: Hog 
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19981 ): L100 lies within Ophiuchus at 125 ± 25 pc fide Geus et al. Ill989h . The results are 
listed in Table [5j If the upper limit kinetic temperature was obtained through the LVG 
analysis, we list it again. 

The derived masses range from 4 to 255 M , depending on the size of the core. The highest- 
mass core, DC267.4-7.5, is also the largest. The smallest and weakest core, DC302.6-15.9, 
does have the lowest mass. A correlation between temperature and mass cannot be estab- 
lished. 

We have found three cores whose masses have been derived from different observations. 



Vilas-Boas et al. I (119941 ) derived the mass of CG30 to be 43 M with an uncertainty of a 
factor of 2.5, but they also quote a shorter distance of 300 pc. With a distance of 450 pc, this 
mass would be 97 M Q , which is within a factor of two of our value (58 M ). iBourke et al 



(119971 ) have calculated the mass of DC297.7-2.8 to be 40 M , using a larger radius and larger 
distance for the source. With the distance we use, 15 pc, this value become s 23 M , which 



is also within a factor of two of our value (12 M ). iKeto fc Myers I (119861 ) calculated the 
mass of DC259.5-16.4 to be 41 M , but they used the much smaller distance of 100 pc. 
Combining their observations with a distance of 450 pc leads to a mass for DC259.5-16.4 of 
830 M . 



4. SUMMARY AND CONCLUSIONS 



In order to study the physical conditions in small isolated molecular clouds, we have mapped 
40 cores in the southern hemisphere in 13 CO(2— >1), 12 CO(4— >3), and 12 CO(7— >6) with 
AST/RO. The sources cover all cores vis ible from the South Pole (5 < —16°) from the 



initial c2d source list (jEvans et al. II2003I ). The AST/RO survey at millimeter and sub- 



millimeter wavelengths provides a direct diagnostic of the dominant gas components in the 
clouds. By observing the 12 CO and 13 CO lines, we have obtained information about the 
physical conditions within the cloud. 13 CO is a tracer of moderately dense (~ 10 3 cm -3 ) 
gas. It is generally optically thin and hence is a good tracer of column density. It probes 
the outer regions of dense cores (10 to 20 K), but not the densest interiors where a small 
percentage of the total 13 CO ma ss is depleted within a small percentage of the total size 
(~ 10 3 AU; iTafalla et al. I (120021 ) ) . The most important submillimeter measurement is the 
12 CO(4— >3) transition: this line is the lowest-lying of the mid- J CO lines that constrain den- 
sity rather strongly. These data have been used to determine meaningful conclusions about 
the physical properties of the molecular cores. The detection of the 12 CO(7— >6) transition 
was unlikely, since the molecular cores are rather cold. From our data, we were able to put 
upper limits on the intensities in this transition, which were important for a Large Velocity 
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Gradient (LVG) analysis. 

We find typical 13 CO(2^1) linewidths to be 1.8 to 2.3 km s _1 for cores with embedded stars 
and 1.6 to 1.9 km s _1 for all others. The slightly higher values for sources with embedded stars 
result from warmer gas whose emission is more subject to thermal and turbulent broadening 
of the lines. Only 10 out of 40 cores have linewidths outside the 1.5 to 2.1 km s _1 window, 
the smallest being 0.8 km s -1 and the largest being 2.7 km s -1 . The largest linewidth is 
found in DC275.9+1.9, a core with an embedded star; the smallest is found in the rather 
weak source DC253.8-10.9, with no signs of star formation. 

The distribution of the 12 CO(4— >3) linewidth is much broader. Typical values lie in the range 
of 1.6 to 2.3 km s _1 for cores with no embedded sources, and 2.6 to 3.7 km s _1 for cores with 
embedded stars. The distribution of their linewidths is rather flat, with the largest linewidth 
of 5.8 km s _1 occurring in DC303.8-14.2. From the 13 CO line intensities, we were able to 
derive the 13 CO column density and the H 2 column density (§ 3). The mean of the 13 CO 
column density is found to be 1.9 x 10 15 cm -2 for cores with embedded stars, and 1.5 x 
10 15 cm~ 2 for all others. Due to its exceptionally strong 13 CO(2^1) line intensity, B59 also 
stands out here with a 13 CO column density of 4.3 x 10 15 cm -2 and an H 2 column density 
of 7.3 x 10 21 cm" 2 . 

For nine cores, we determined the kinetic temperature and molecular hydrogen density 
through an LVG analysis. For the 12 CO(7— >6) transition, we determined upper intensity 
limits, since we have no significant detection of this line. As a result, all temperatures and 
densities are also upper limits. In addition, we give a lower temperature limit, assuming the 
absence of 12 CO(7— >6) emission. For these nine cores, four of which have embedded stars, 
we find an average lower temperature limit of 16 K and an average upper limit of 26 K. The 
lower temperature limit is smaller than 21 K for all cores. The average upper limit on the 
molecular hydrogen density was found to be 10 3 3 cm -3 for all cores. The average lower limit 
on the molecular hydrogen density was found to be 10 2 ' 9 cm -3 . Due to the small sample 
size, we cannot detect a significant difference between the mean densities of cores with and 
without embedded stars. 

We have obtained masses for nine cores. They range from 4 to 255 M , depending on the size 
of the core. A correlation between temperature or density and mass cannot be established, 
which is mainly due to the small sample size. 

In general, the sources in this survey behave rather uniformly with only very few exceptions. 
Cores with embedded sources tend to have larger linewidths and higher 13 CO column density 
than cores with no embedded stars. 



From the comparison to previously published CO surveys of molecular cores (IKeto fc Myers 
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19861 ; IVilas-Boas et al. Ill994j . 12000). we can say that the c2d sample of molecular cores has a 
range of properties (linewidth, column density, size, mass, embedded stars) similar to those 
of past studies. Only three sources from the sample covered in this survey have previously 
been observed in other CO surveys, albeit in lower rotational transitions. To obtain even 
more accurate results with an LVG analysis, a survey of these sources at one other CO 
transition would be very valuable. 

The authors would like to thank Jacob Kooi of Caltech for his excellent work on the telescope 
and the receivers. We thank the receiver group at the University of Arizona and Greg Wright 
of Antiope Associates for their work on the instrumentation; R. Schieder, J. Stutzki, and 
colleagues at the University of Koln for their AOS's. We are grateful to S. Kim and M. Yan 
for making their LVG code available to us. We acknowledge the use of NASA's SkyView 
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Table 1. Source list. 
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Table 1 — Continued 



Source name 


Other names 
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5 [J2000] 1 
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1 Coordinates of the map center 

2 Cores with known embedded sources as detected using IRAS are marked (*). 

3 12 CO(4^3) data on DC297.7-2.8 will be presented by Bourke et al., in preparation 



Table 2. Observed and derived parameters of the 13 CO(2^1) lines. 
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Table 2 — Continued 
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1 Coordinates of spectrum with highest integrated intensity. 

2 Linewidth obtained using high resolution spectrometer data. 
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Table 3. Observed parameters of the 12 CO(4— >3) lines. 
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n 


A 


0.4 


DC266.0-7.5 


08 


21 


13.6 


-49 


54 


47 


3.8 


2 

Zj . VJ 


± 


0. 


.2 


4.9± 


n 

\J ■ 




0.3 


* DC267.4-7.5 


08 


25 


45.1 


-50 


59 


18 


5.5 


1 8 


± 


0. 


1 


10.9 


± 


n 

VJ 


,2 


0.2 


DC267.5-7.4 


08 


26 


56.1 


-51 


04 


31 


6.0 


2 7 


± 


0. 


.2 


5.9 


± 


n 

VJ 


■ • ) 


0.3 


DC253.6+2.9 


08 


28 


58.4 


-33 


47 


12 


6.3 


1 3 

X . *J 




0. 


.1 


6.4 


± 


n 

VJ 


3 


0.3 


DC267.6-6.4 


08 


32 


03.9 


-50 


32 


21 


5.8 


1 6 

X . V ' 


± 


0, 


1 


5.2 


± 


n 

VJ 


3 


0.3 


DC269.4+3.0 


09 


22 


16.3 


-45 


47 


25 


-1.4 


2 

_ . VJ 




0. 


.2 


3.3 


± 


n 

VJ 


2 


0.3 


DC274.2-0.4 


09 


28 


50.2 


-51 


36 


40 


5.9 


2 2 


± 


0, 


.2 


2.6 


± 


n 

VJ 




0.2 


DC272.5+2.0 


09 


31 


11.9 


-48 


37 


52 


-3.8 


2 7 


± 


0. 


.9 


1.5 


± 


n 

VJ 


2 


0.3 


* DC275.9+1.9 


09 


46 


30.1 


-51 


01 


12 


-5.7 


2 6 


± 


0. 


3 


4.2 


± 


n 

VJ 


,4 


0.4 


DC291.0-3.5 


10 


59 


02.4 


-63 


40 


43 


-4.8 


1 

X . W 


± 


0. 


.1 


5.0 


± 


n 

VJ 


,4 


0.6 


DC291. 1-1.7 


11 


06 


13.6 


-62 


03 


13 


-4.3 


1.2 


± 


0. 


.1 


3.5 


± 


n 

VJ 


,2 


0.3 


DC300.2-3.5 


12 


24 


19.0 


-66 


12 


52 


-5.7 


i i 
i . i 


it 


0. 


1 


2.7 


± 


0. 


.2 


0.2 


Mu8 


12 


28 


07.6 


-71 


19 


05 


3.25 


2.5 


± 


0. 


3 


6.4 


± 


0. 


.5 


0.7 


DC300. 7-1.0 


12 


31 


29.5 


-63 


44 


32 


-4.9 


3.5 


± 


0, 


1 


5.9 


± 


0. 


,1 


1.3 


DC301.2-0.4 


12 


36 


19.4 


-63 


11 


09 


-4.7 


2.6 


± 


0. 


.5 


2.4 


± 





.3 


0.3 


DC302.6-15.9 


12 


45 


38.7 


-78 


47 


21 


3.8 


1.7 


± 


0. 


1 


2.2 


± 


0. 


.1 


0.5 


* DC302.1+7.4 


12 


45 


38.5 


-55 


24 


48 


-14.7 


2.2 


± 


0. 


.3 


3.7 


± 





,4 


0.4 


* DC303.8-14.2 


13 


07 


26.8 


-77 


01 


35 


3.4 


5.8 


± 


0. 


.5 


6.7 


± 





.3 


0.2 


DC327.2+1.8 


15 


42 


32.4 


-52 


48 


09 


-0.8 


3.2 


± 


0. 


,5 


3.1 


± 


0. 


.3 


0.3 


DC339.1 + 11.7 


15 


59 


31.4 


-37 


36 


57 


5.1 


1.3 


± 


0. 


1 


6.6 


± 


0. 


.5 


0.5 


DC338.6+9.5 


16 


04 


33.2 


-39 


38 


56 


4.4 


1.0 


± 


0. 


1 


3.6 


± 


0. 


.3 


0.3 


DC337.6+7.6 


16 


07 46.2 


-41 


43 


49 


4.7 


1.7 


± 


0. 


1 


9.5 


± 





.1 


0.5 


* L43 


16 


34 


25.8 


-15 


48 


27 


0.4 


3.0 


± 


0, 


3 


14.8 


± 


1 


.0 


1.0 


DC346.0+7.8 


16 


36 


57.2 


-35 


36 


38 


3.3 


1.8 


± 


0. 


.2 


6.8 


± 


0.4 


0.4 
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Table 3 — Continued 



Source a [J2000] 1 5 [J2000] 1 v LSR Av J T MB (v) dv rms Noise 





pi m 


1 


[° 


'"] 


[ km s x ] 


[ km s 


- 1 ] 


[K km s 


- 1 ] 


[K] 


DC346.3+7.8 


16 


37 


47.4 


-35 


28 


53 


3.4 


1.8 


± 


0.1 


7.1 


± 


0.3 


0.3 


DC346.4+7.9 


16 


37 


57.5 


-35 


15 


11 


3.3 


1.6 


± 


0.1 


4.8 


± 


0.2 


0.3 


L63 


16 


50 


04.7 


-18 


01 


04 


5.8 


1.2 


± 


0.1 


6.7 


± 


0.2 


1.4 


* CB68 


16 


57 


03.8 


-16 


03 


32 


5.6 


3.4 


± 


0.4 


7.5 


± 


0.4 


1.2 


* B59 


17 


11 


17.6 


-27 


26 


28 


4.0 


3.6 


± 


0.3 


18.1 


± 


0.7 


0.9 


★ L100 


17 


16 


16.0 


-21 


01 


51 


1.4 


1.8 


± 


0.3 


14.2 


± 


1.6 


0.4 


LI 772 


17 


19 


50.6 


-26 


44 


40 


3.5 


2.1 


± 


0.2 


5.1 


± 


0.3 


0.4 


B72 


17 


23 


04.7 


-23 


42 


28 


5.1 


1.9 


± 


0.3 


6.0 


± 


0.5 


0.8 


L328 


18 


17 


02.1 


-18 


04 


24 


6.9 


1.2 


± 


0.3 


3.3 


± 


0.7 


0.9 



1 Coordinates of spectrum with highest integrated intensity. 

2 Cores with known embedded sources as detected using IRAS are marked (*) . 



Table 4. Results from the LVG model. 



Source 


a [J2000] 1 


5 [J2000] 1 


nil 


rpl2 


rpl2 


T l2 


ni 3 

1 2->1 


T 


T 

- 1 - kin 


log n(H 2 ) 


log n(H 2 ) 




rh m si 
V J 


L J 


L J 


[Kl 

L J 


L J 






[Kl 

L J 


[Kl 

L J 


[cm -3 l 

L J 


fcm _3 l 

L J 












u.l. 3 


U.l. 3 




u.l. 3 


l.l. 3 


u.l. 3 


l.l. 3 


DC259.5-16.4 


07 14 03.2 


-48 27 40 


1.6 


0.8 


0.5 


0.6 


0.5 


22 


14 


3.2 


2.9 


CG31A-C 


08 08 49.2 


-35 57 51 


3.1 


4.3 


1.4 


0.3 


1.4 


28 


17 


3.2 


2.8 


* 2 CG30 


08 09 26.3 


-36 03 51 


1.6 


1.8 


0.8 


0.4 


1.2 


30 


18 


3.3 


2.9 


★ DC267.4-7.5 


08 25 48.3 


-50 59 48 


2.5 


2.8 


0.3 


0.1 


1.1 


22 


14 


3.2 


2.9 


DC267.5-7.4 


08 26 46.5 


-51 04 01 


1.6 


1.5 


0.4 


0.3 


0.9 


26 


15 


3.2 


2.9 


DC274.2-0.4 


09 28 47.0 


-51 35 10 


1.0 


0.6 


0.2 


0.4 


0.6 


18 


12 


3.7 


3.0 


DC339.1+11.7 


15 59 11.2 


-37 36 57 


2.3 


3.4 


1.3 


0.4 


1.5 


34 


19 


3.1 


2.8 


* B59 


17 11 19.8 


-27 25 58 


4.0 


4.1 


1.9 


0.5 


1.0 


28 


16 


3.3 


2.9 


* L100 


17 16 07.4 


-20 57 21 


1.7 


2.3 


0.6 


0.3 


1.4 


30 


20 


3.1 


2.7 



1 Coordinates of position where LVG analysis was performed. 

2 Cores with known embedded sources as detected using IRAS are marked (*). 

3 Upper and lower limits are denoted u.l. and 1.1. 
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Table 5. Masses derived from the H2 column densities. 



Source 


Distance 1 


Radiusi 2 


Radius2 2 


N(H 2 ) 




Mass 






[pc] 


[pc] 


[P C ] 


[10 21 cm 




[M ] 


u.l. 3 [K] 


DC259.5-16.4 


450 ± 50 


0.9 


0.7 




1.3 


48 ± 12 


22 


* 4 CG30 


450 ± 50 


0.8 


0.7 




1.8 


58 ± 14 


30 


* DC267.4-7.5 


450 ± 50 


1.4 


1.2 




2.9 


255 ± 60 


22 


DC253.6+2.0 


450 ± 50 


0.9 


0.7 




2.4 


103 ± 26 




DC274.2-0.4 


450 ± 50 


0.8 


0.7 




1.5 


56 ± 14 


18 


DC297.7-2.8 


150 ± 30 


0.2 


0.3 




3.3 


12 ± 4 




DC302.6-15.9 


150 ± 30 


0.2 


0.3 




0.9 


4 ± 2 




* DC303.8-14.2 


178 ± 18 


0.5 


0.4 




1.5 


19 ± 4 




* L100 


125 ± 25 


0.4 


0.4 




1.9 


14 ± 6 


30 



Woermann et al. 


f 


2001 


I); 


Franco 


f 


1995 


); 


Corradi et al. 


Knude & Hop 


f 


1998): 


ie Geus et al. ( 


198S 


) 



2 Radii of the elliptical half power width of the core. 

3 Upper limit as derived from the LVG analysis. 

4 Cores with known embedded sources as detected using IRAS are marked (*). 
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Fig. la. — 13 CO(2^1) (left) and 12 CO(4— >3) (right) spectra with the highest integrated 
intensities. Coordinates of the spectra are given in Tables [2] and [3J 
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Fig. 2.— Histograms of the 13 CO(2^1) linewidth (left) and 12 CO(4^3) (right). Grey bars 
represent sources with no embedded stars, black bars represent sources with embedded IRAS 
sources. 
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Fig. 3. — Histogram of the 13 CO column density for all sources. Grey bars represent sources 
with no embedded stars, black bars represent sources with embedded IRAS sources. 
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Fig. 4a. — Maps in 2 transitions of sources chosen for LVG analysis: 13 CO(2^1) (left/top) 
and 12 CO(4— >3) (right/bottom). Contours are in steps of 10% of the peak intensity, listed in 
each map. The effective beamsizes are indicated by the circles and are 216" in 13 CO(2^1) 
and 117" in 12 CO(4^3). (a) DC259.5-16.4, DC274.2-0.4, DC339.1+11.7, (b) L100, (c) 
CG31A-C, CG30, (d) DC267.5-7.4, DC267.4-7.5, (e) B59 
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Fig. 5.— Maps of DC253.6+2.0, DC297.7-2.8, DC302.6-15.9, and DC303.8-14.2 in 
13 CO(2^1). The effective beamsize of 216" is indicated by the circle. 
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